Baryon Acoustic Oscillation Intensity Mapping as a Test of Dark Energy 
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The expansion of the universe appears to be accelerating, and the mysterious anti-gravity agent 
of this acceleration has been called "dark energy" . To measure the dynamics of dark energy, Baryon 
Acoustic Oscillations (BAO) can be used. Previous discussions of the BAO dark energy test have 
focused on direct measurements of redshifts of as many as 10 9 individual galaxies, by observing the 
21cm line or by detecting optical emission. Here we show how the study of acoustic oscillation in the 
21 cm brightness can be accomplished by economical three dimensional intensity mapping. If our 
estimates gain acceptance they may be the starting point for a new class of dark energy experiments 
dedicated to large angular scale mapping of the radio sky, shedding light on dark energy. 



Introduction. — To understand dark energy and 
sharply test theories of its character, it is necessary to 
precisely measure the last half of the expansion history. 
This period corresponds to the redshift range < z < 2. 
Many techniques have been proposed to study the late- 
stage expansion history, and some of the most promising 
make use of Baryon Acoustic Oscillations (BAO)[l|. In 
addition to producing CMB structure, the acoustic oscil- 
lations also produced density structure in the atomic gas 
and dark matter, which is still detectable today. Many 
groups have reported 2 to 3 a detections, in the low red- 
shift universe, of periodic structures in the density of 
galaxies at the predicted wavelengths across the skyQ. 
Because the acoustic waves are frozen in after recom- 
bination, the BAO peak wavelengths can be used as a 
cosmological standard ruler: observation of the angular 
size of the peak wavelengths across a range of redshifts 
allows accurate measurement of the expansion history. 

In this letter we present calculations of structure in 
the three dimensional brightness due to the hypcrfinc 
transition of neutral hydrogen (HI) at 21 cm wavelength. 
We show that, via 21 cm emission, baryon oscillations 
could be precisely measured, using a telescope just 200 
wavelengths in diameter, since each cosmic cell of the 
appropriate scale contains more than 10 12 M Q of emit- 
ting neutral hydrogen gas. We present forecasts on the 
dark energy constraints from this new type of observa- 
tion, which are competitive with the best proposed dark 
energy experiments [3j. Where needed in the analysis we 
adopt WMAP3 values for the cosmological parameters[4| , 
n m = 0.24, fl b = 0.04, n A = 0.76, and h = 0.73, where 
fi m , fib, and Oa are the matter, baryon, and dark energy 
fractional density, respectively, and h is the dimension- 
less Hubble parameter. We also use the WMAP3 error 
limits for these values. 

A recent paper proposed a test of dark energy mod- 
els by measurement of baryon oscillations in the 21 cm 
brightness field at redshifts z > 3 |. Using only 
such high redshift data two cosmological parameters — 
departures of spatial flatness — 1 and a slow change of 
dark energy equation of state — are nearly indistinguish- 



able. It is only when precise flatness £1 = 1 is assumed 
that such high redshift data can be used to constrain dark 
energy models. In this letter we use the more standard 
approach and assume that ft is not perfectly measured. 
Then data at lower redshift become essential if dark en- 
ergy models are to be tightly constrained. 

Detectability. — To measure the oscillations it is not 
necessary to detect individual galaxies, but only to mea- 
sure the variation in HI mass on large scales. A familiar 
analogy is the study of galaxies via images of their sur- 
face brightness. Except for the nearest of galaxies, we do 
not detect the individual stars that produce the optical 
emission, indeed the pixels of most galaxy images repre- 
sent the emission of millions of stars. Never-the-less such 
images allow study of galactic structure. 

Above the third BAO peak nonlinear evolution atten- 
uates BAO structure, so the third peak is at the smallest 
spatial scale we need consider. This peak has a wave- 
length 35ft. -1 Mpc. A Nyquist sampled map therefore 
needs pixels of size 18/i _1 Mpc. At z — 1.5 the corre- 
sponding angular wavelength is 20 arc minutes, which 
requires a telescope of approximately 200 wavelength, or 
100 meters, to resolve. 

We base our estimate of neutral hydrogen mass con- 
tained in a 18ft, -1 Mpc volume ( M H i ~2x 10 12 M Q ) on 
Lya absorption studies of the cosmic density of neutral 
gas [6|, which find Ohi ~ 1 X 10~ 3 at z ~ 1, with a weak 
dependence on redshift. 

The average sky brightness temperature due to the 
21cm line is about 300/iK. The mean brightness tem- 
perature Tf, of the 21cm line can be estimated usingQ 




Most observation of the sky are measurements of varia- 
tions in temperature AT], ~ T{,5, where 1 + 8 = p g / p g is 
the normalized neutral gas density, and et=(l + z) -1 is 
the scale factor. 



On I8h 1 Mpc scales spatial variations of the sky 
brightness due to 21 cm emission in the cosmic web 
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FIG. 1: The power spectrum of neutral hydrogen hyperfine 
emission at z = 1.5 is shown by the solid line. The dashed line 
indicates the map noise from a 200-by-200-meter radio cylin- 
der telescope operated in transit mode for 100 days, assuming 
a telescope system temperature of 50 K. The horizontal dot- 
ted line indicates the non- linear scale at z = 1.5. The vertical 
lines indicate the location of peaks in Figure [2] 



FIG. 2: Baryon Acoustic Oscillations. To show BAO we plot 
the ratio of matter power spectrum that contains contribu- 
tions from baryons, calculated using a fitting formula 1 101]. over 
a power spectrum with no baryons. To single out the oscil- 
lation, we further divide the ratio by a smooth fitting curve 
to the overall spectrum. The error bars represent projections 
of the sensitivity possible with a 200-by-200-meter cylinder 
telescope. 



structure presents a sky noise at the 150 /xK level. Fig- 
ure[T]shows the amplitude of this cosmological large scale 
structure, while Figure [5] plots the small oscillations im- 
printed by BAO. The variation of gas density due to 
baryon oscillations is small compared to the fluctuation 
of the large-scale structure. To detect the BAO signals 
beneath this noise Fourier analysis of a large regions of 
sky is required. The observations should cover as much 
sky as possible. 

For a filled aperture telescope the RMS variation in 
measured sky temperature AT = T sys /y/ Af t- lnt , where 
AT is the temperature in Kelvin, T sys is the system 
temperature of the radio telescope in Kelvin, Af is the 
observing bandwidth in Hertz and tint the integration 
time in seconds. Assuming current technology, which 
allows T sys = 50K, one needs an observing bandwidth- 
time product of 2 x 10 11 to bring the thermal map noise 
down to ~ 100/iK, below the cosmic structure noise. In 
frequency space 18 Mpc spans roughly 3 MHz at 
z = 1.5, so the required observation time is ~ 6 x 10 4 
seconds, about 18 hours per pixel. Small regions of sky 
might be mapped at this noise level using existing tele- 
scopes such as the Green Bank Telescope (GBT), the 
Effelsberg 100-meter Radio Telescope or the Parkes Tele- 
scope. 

Brightness sensitivity is not improved by increase of 



telescope aperture, instead multiple receivers, either on 
a single reflector or on individual reflectors, can be used 
to increase the accumulated integration time. 

For sensitivity estimates, we consider an intensity map- 
ping telescope, dubbed IM, covering a square aperture of 
size 200m x 200m. We subdivide this area into 16 cylin- 
drical sections, each 12.5m wide and 200m long, analo- 
gous to the Hubble Sphere Hydrogen Survey (HSHS)(llj. 
The cylinders are aligned North-South allowing the tele- 
scope to instantaneously sample almost the entire merid- 
ian. This geometry allows all-sky coverage each 24 hour 
period in a telescope with no moving parts. The esti- 
mated sensitivity of the survey is shown in Figure [T] 

Foregrounds. — The desired BAO signal is smaller than 
other sources of sky brightness structure. The largest 
interference is due to synchrotron emission from the 
Galaxy. The spatial structure of this emission exceeds 
the BAO signal by 10 3 . Similar emission by extragalac- 
tic radio sources distributed across the sky will contribute 
sky brightness fluctuations on smaller scales. However 
all these foregrounds are spectrally smooth, with bright- 
ness temperature T oc f a , so one expects to be able to 
subtract them spectrally. The image from each 3 MHz 
slice of the frequency spectrum contains a unique BAO 
pattern, largely uncorrelated with the neighboring slices. 
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In contrast, the foregrounds and backgrounds are highly 
correlated. By subtracting the average of both neigh- 
boring slices the BAO signal can be retained while the 
other emission is removed. Of course this subtractions 
will be imperfect, because of imperfections in the tele- 
scope and spatial variations in the astrophysical source 
spectrum. A similar challenge is faced by the higher- 
redshift 21cm experiments probing the Cosmic Reioniza- 
tion era and many foreground subtraction methods have 
been proposed ([l2j and references therein). 

The spectral index a for Galactic emission varies across 
the sky. The frequency-double-difference technique de- 
scribed above suppresses the Galactic sky brightness but 
leaves a residual of order AT = Aa(Af / f) 2 T to first 
power in Aa, where Aa is the uncertainty in spectral in- 
dex and A/ the frequency difference. Using the WMAP3 
measurement we estimate 
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where I is the spherical harmonic number. The normal- 
ization factor of 1.3 Kelvin temperature fluctuation is de- 
rived by extrapolating from the WMAP3 measurement of 
Galactic foreground variation at 20 GHz on 3.7° degree 
scale, along with the spectral index of synchrotron emis- 
sion at mid-Galactic latitude 13]. Adopting Aa = 1, a 
conservative upper limit, we estimate a residual Galactic 
interference signal which is plotted in Figure [3] 

Figure [3] shows various potential limitations to BAO 
measurements using 21 cm sky brightness mapping. The 
top exclusion comes from the non-linear scale, which 
erases the baryon wiggles. The left exclusion comes from 
the finite volume, which sets a lower bound on the low- 
est wavenumber that fits within a given volume. The 
lower right is a rough estimate of the impact of fore- 
grounds using the double frequency difference described 
above. The shaded region reflects the region where the 
estimated residual foreground from Eqn. exceeds the 
expected signal in Eqn. |T]). 

Parameter Estimation. — To estimate the sensitivity of 
this technique for measurement of dark energy parame- 
ters, we follow the procedure in Seo and Eisenstein (2007) 
14] and carry out the standard figure-of-merit calcula- 
tion used in the DETF report 0. Our calculations in- 
clude the effects of redshift space distortion and a sim- 
plified model of non-linearity using a smoothing window. 
This type of analysis, the conversion from survey data 
to cosmological parameters, is an area of active research, 
and several effects such as mode-mode coupling are still 
being investigated. Our use of shells of 21 cm surface 
brightness introduces two new types of unknowns into the 
analysis: the HI density, and the bias of 21 cm sources. 
To accommodate these uncertainties we make one sensi- 
tivity estimate using conservative values for the unknown 
quantities (IM) as our baseline model, and one with a 



0.1 



o 
a 
S 

\ 



0.01 - 




1 1.5 
redshift 



FIG. 3: The observable parameter space in redshift and in 
scale (k) for BAO. The shaded regions are observationally 
inaccessible (see text). The horizontal lines indicate the scale 
of the first three BAO wiggles, and the dashed lines show 
contours of constant spherical harmonic order i. 



worst case assumption that each of these effects works to 
increase the error (IMW) . The IM case estimate assumes 
that l^Hi remains at the present day value of 0.0005, in- 
dependent of redshift, and we adopt a bias 6=1. The 
shot noise from present day HI galaxies [HI is negligible. 
Motivated by the reconstruction of BAO degradation due 
to non-linearity in [l6| the Seo and Eisenstein damping 
scale is shortened by a factor of two. For the worst case 
(IMW) estimate we let fl H i to be 0.0004. We adopt a 
bias b = 0.8, which is the current estimate of the low 
z value. This increases the noise, and seems securely 
to be a worst case because galaxies are unlikely to be 
both antibiased and rare. We assume that 21 cm sources 
are rare and massive at high redshift, and we use the 
Lyman-Break Galaxy luminosity function (with a den- 
sity of 0.004 h 3 Mpc -3 ) for the shot noise calculation. 
The Seo and Eisenstein damping scale is assumed. 

The time variable dark energy density is parametrized 
as pde oc (1 + z) 3 ( 1+w ^; an ordinary cosmological con- 
stant corresponds to w — — 1, and the principal goal of 
dark energy experiments is to measure deviations of w 
from — 1. We first determine the Fisher matrix for the 
tangential and radial distances errors. From these, we 
compute the errors on dark energy parameters. The re- 
sulting projected errors to the equation of state are shown 
in Figure |4] Table H] quantifies the errors, which can be 
compared to other dark energy surveys. 

The mapping of cosmological large scale structure by 
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TABLE I: Fiducial cosmological parameters and error bars for combining the Planck mission with the fiducial BAO intensity 
mapping experiment such as IM (P+IM), with the Dark Energy Task Force 0] (DETF) Stage II projects (P+II), plus the Stage 
III projects (P+II+III), plus the Stage IV projects (P+II+III+IV), and all above combined (P+II+III+IV+IM). Here IM is 
the conservative best guess for intensity mapping experiment parameters, while IMW is the worse case scenario (see text). 
DETF Stage II is well-defined which includes weak lensing (WL), cluster evolution (CL), and supernova (SN) experiments; 
for Stage III we use the optimistic photometric versions of WL, CL, and SN; for Stage IV we use the optimistic space-based 
WL, CL, and SN experiments. The table entries are shown using the units of the DETF Report: IM compares favorably with 
other proposed Stage III projects. The dark energy equation of state is parametrized as w — w p + (a p — a)w' . Qoe and flfc 
are the dark energy and curvature density, respectively. n s is the spectral index of the primordial power spectrum, and A its 
normalization. The Figure of Merit (FoM) refers to the inverse volume in the 95% confidence ellipse. 
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FIG. 4: The 1-er contour for IM combined with Planck (inner 
thick solid for baseline model, outer thin solid for worst case), 
the Dark Energy Task Force Stage II projects with Planck 
(outer dotted), the Stage II and III projects with Planck (in- 
termediate dotted), the Stage II, III and IV projects with 
Planck (inner dotted), and all above experiments combined 
(dashed, again thick for baseline, thin for worst case; the two 
contours are nearly indistinguishable). 



three dimensional 21 cm intensity mapping is an ap- 
proach that can be tested with existing telescopes today, 
which appears to offer very competitive constraints on 



dark energy parameters. This new tool may soon allow 
astronomers to economically map huge volumes of the 



universe. 
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